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The invariant form of the hyperfine interaction between muhipolar moments and the 
nuclear spin is derived, and apphed to discuss possibihties to identify the antiferro-octupolar 
(AFO) moments by NMR experiments. The ordered phase of Np02 and the phase IV of 
Cei_a;Laa;B6 are studied in detail. Recent ^^O NMR for polycrystalline samples of Np02 are 
discussed theoretically from our formulation. The observed feature of the splitting of ^''O 
NMR spectrum into a sharp line and a broad line, their intensity ratio, and the magnetic field 
dependence of the shift and of the width can be consistently explained on the basis of the 
triple q AFO ordering model proposed by Paixao et. al. Thus, the present theory shows that 
the ^^O NMR spectrum gives a strong support to the model. The 4 O sites in the fee Np02 
become inequivalent due to the secondary triple q ordering of AF-quadrupoles: one cubic 
and three non-cubic sites. It turns out that the hyperfine field due to the antiferro-dipole 
and AFO moments induced by the magnetic field, and the quadrupolar field at non-cubic 
sites are key ingredients to understand the observed spectrum. The controversial problem 
of the nature of phase IV in Cei-ajLa^jBg is also studied. It is pointed out that there is a 
unique feature in the NMR spectra, if the ^^{T^ = T^ = T^) AFO ordering is realized in 
Cci-xLaa-Bg. Namely, the hyperfine splitting of a B atom pair on the (i, ^, ±u) sites crosses 
zero on the (110) plane when the magnetic field is rotated around the [001] axis. 

KEYWORDS: antiferro-octupolar order, octupolar moment, hyperfine interaction, NMR, 
Np02, CeBg, antiferro-quadrupolar order 



1. Introduction 

The antiferro-octupolar (AFO) moments, which are odd in time-reversal symmetry but 
are distinct from dipole moments, are now widely recognized as a "hidden" order parameter.^ 
The AFO order is a promising candidate to resolve the breaking of the time reversal symmetry 
without the antiferromagnetic dipole (AFM) order. However, a direct observations of the AFO 
is very difficult. 

There are two ways to realize the AFO moments. One is induced AFO moments in the 
antiferro-quadrupolar (AFQ) state by applying the static magnetic field. ^ Another is a spon- 
taneous AFO ordering, which has been recently proposed for Np02 and Cei-j^LaxBg.^"^ In 
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this paper, we wish to study a possibihty to detect the AFO moments in NMR experiments 
through the hyperfine interaction. For this purpose we derive a general form of the hyper- 
fine interaction between the multipolar moments and the nuclear spin. Actually we studied 
this problem before for CeBg. The approach was to derive a phenomenological form of the 
hyperfine interaction from the symmetry.^' ^'^ However, there is a simpler and equivalent ap- 
proach, which we use in this work. It is a straightforward generalization of rcf.ll , in which 
the invariant coupling form between multipolar moments was derived by using a method of 
symmetrized molecular orbital theory. The usefulness of this approach to discuss the hyperfine 
interaction will be demonstrated in this paper. 

To have octupolar moments as independent local degrees of freedom, the material must 
have high symmetry such as cubic. There are three cubic systems for which the AFO order 
of some sort is either realized or is likely to be realized. The first system is CeBg, which 
has been studied extensively.^*^' ""^^"^"^ There are a number of anomalous features in the phase 
diagram as well as in the nature of the low-temperature order phases, phase II in particular. 
They are naturally explained by noting the interaction effect between the field-induced AFO 
moments, especially of the Txyz type type).^*^' ""^^ However, a crucial progress was achieved, 
when the apparent inconsistency between the neutron diffraction-*^^ and NMR^'' for the phase 
II was successfully resolved by recognizing that the hyperfine coupling is present between the 
B nuclear spin and the induced T^yz AFO moment of Ce ion.^'^*^'^^ This is the first case in 
which the contribution from AFO moments to the hyperfine coupling has been discovered. 
Therefore the NMR is now regarded as a direct proof for the field-induced AFO moment in 
CeBg. 

The second system is Np02; the nature of the order in Np02 below 26K was a mystery for 
a long time. Carrying out a resonant X-ray scattering, Paixao et. al. have recently suggested 
that their results can be explained well by assuming the triple q AFQ order of Fs type.^ They 
proposed the triple q AFO as the primary order parameter because any dipole order have 
not been observed though Np''^ is the Kramers ion. The triple q AFQ is induced from the 
primary triple q AFO. A new development on this problem has been brought by Tokunaga et. 
aVs ^^O (/ = 5/2) NMR experiment for polycrystalline samples. They found that the O ion 
sites, which are equivalent in the high temperature phase (Fm3m), split into two inequivalent 
groups. The experimental results have the following unique features: 

(I) 1/4 of the^^O nuclei contributes to a sharp resonance line, while the rest 3/4 gives a broad 
resonance line. 

(II) The shift of the sharp line is approximately proportional to the strength of the applied 
magnetic field. 

(III) The width of the broad line is given as the sum of a constant and a part proportional to 
the strength of the magnetic field. 
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(IV) The shift of the sharp hne and the magnetic field induced width of the broad Hne have 
comparable magnitude. 

We show in this paper that those experimental results can be naturally explained by the 
quadrupolar field and the hyperfinc field splitting in the triple q state. The hyperfine field is 
caused by the AFO and/or AFM moment induced by the application of the magnetic field. 
Therefore the ^'^O NMR is a fingerprint of the triple q order of AFO and AFQ. 

The third system is Cex^^i-x^&- It is known that different from pure CeBg a new phase, 
which is called phase IV, appears instead of phase II in the low magnetic field region of 
Ce^^Lai^ajBg {x ~ CTS).^*^"^^ The nature of phase IV remains controversial. The time reversal 
symmetry breaking is reported for this phase. ^^'^^ However, no indication of the AFM or- 
der has been detected experimentally, although neutron diffraction studies were carried out 
extensively.^^ A characteristic cusp of the magnetic susceptibility has been observed at the 
transition temperature between phase I and phase IV. From these results together with other 
experimental indications,^''"^^'^^'^'' an AFO order of the Fs type has been recently proposed 
for the phase IV. A recent NMR experiment indicated a broadening of NMR lines in the 
phase IV, showing a sign of time reversal symmetry breaking. ^"^ No further information is 
obtained from this NMR experiment. Quite recently, the importance of the uniaxial stress 
to remove the ambiguities due to the domain structure has been demonstrated. Morie et. al. 
pointed out that the magnetic susceptibility shows an anisotropy in contrast to the isotropic 
behaviors reported previously for ambient pressure. This result is qualitatively similar to 
what one obtains theoretically from the AFO of the F5 type with [Tx) = (Ty) = (Tz)- Here, 
Tx , Ty and Tz are the three components of the Fs type octupolar moment. The NMR study 
under uniaxial stress is desirable to clarify the situation for phase IV. In this paper we will 
point out a possibility to identify the T^^ type AFO order through the NMR experiment. We 
will show that if the {Tx) = {Ty) = {Tz) type order is realized as Kubo and Kuramoto, and 
Morie et. al. assumed,^' a characteristic feature of the NMR splitting should be present in 
the field-direction dependence of the B NMR from {\,\, ±m) sites. 

This paper is organized as follows. In §2, the hyperfine interaction with multipolar moments 
of 4f electrons is derived for Np02 on the basis of new method. In §3, the recent ^^O NMR 
experiment on Np02 is studied theoretically to show that it is consistent with the triple q 
structure of AFO and AFQ. In §4, the hyperfine interaction in CeBg is re-derived. The field 
direction dependence of the NMR line is discussed in §5, having the phase IV problem in our 
mind. 

2. Invariant Coupling in Np02 and Splitting of O sites into Non-equivalent Sites 

The Np ions in Np02 form the f.c.c. lattice. The position vector p of the eight O ions in 
the f.c.c. cube is given in Table II, where the length of the edge of the cube is chosen as 2a. 
The local symmetry around the O sites has the character. Prom the discussions described 
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in Appendix A, the hyperfine interaetion of the ^''O nuclear spin (/ = 5/2) located at p with 
multipolar moments of Np 5f electrons is given by 

Hhdp) = ^ Ixip) ci^i-^Jxiq){cxCyCz-isg{p)sxSySz) 

+Cl,2-4={^2/(g )(iS5'(p )SzCxCy - CzSxSy) 



16 

+Jz{q ) {}sg{p )syCzCx - CyS^Sx) } 

+ci,3;^{Tf )szCxCy - c^SxSy) 

-TziQ ){^sg{p )syCzCx - CySzSx))^ 
4 -1 

+Ci^i-^xyz{q){}Sg{p )SxCyCz - CxSySz) 

+ (cyclic permutation of x^y and z) . (1) 

Here Qj-'s are coupling constants, and sg{p) = 2)/", and Cu{su) represents 

cos(qi,/2) (sin(Q'iy/2)). The summation over q is assumed in eq. (1). Under the cubic rotation 
group the octupolar operator of r4 type has the same symmetry character with the dipole 
operator. Therefore the dipole operator J in eq. (1) should be read as a linear combination of 
the pure dipole and the octupolar operator of T4 type. The terms containing the factor isg{p ) 
are due to the T^ site symmetry around each O ion. The imaginary factor in this quantity is 
only a seeming, which is compensated by the factor e'^ P . 

Let us assume the AFO order of triple q type: Ti{R) = T^e^Q , T^{R ) = 
tI^q^Q yR and T^{R) = T^'Q , with = 7r(l,0,0), Qy = 7r(0,l,0) and Q ^ = 
7r(0, 0, 1), as assumed in ref. (4 ). Then we easily find from eq. (1) that the AFO moments in 
this order do not induce any hyperfine field on the O atoms. 

As noted in refs. (4) and (11) and seen from Table I, the triple q ordering induces the AFQ 
of the Fg type: Oyz{R ) = Oyze'Q , OzxiR ) = 0,,e'^ and O^yiR ) = Oxye'Q . If 
Oyz = Ozx = Oxy = hold, the cubic symmetry of the crystal is not broken. 

Let us consider the quadrupolar interaction between the ^^O nuclear moment on p and the 
quadrupolar moment of Np ion. The interaction form is given by following the discussion in 
Appendix B: 

Hqqip ) = e'PI C2,i^((0„(p )Ou,q + Oy{p )Oy^q ){cxCyCz - \sg{p )sxSySz)) 

+ C2,2^^{\{-Ou{p ) + V?>Ov{p ))Oyz,q {-'^Sg{p )SxCyCz + CxSySz) 

+Ou{p )Oxy,q {-'iSg{p )SzCxCy + CzSxSy) 

+\{-Ou{p ) - VSOyip ))Ozx,q {-isg{p )syCzCx + CySzSx)) 
+ C2,3^{Oyz{p )l{-Ou,q + V^Oy^q )iisg{p )sxCyCz - CxSySz) 
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+Ozx{p )^i-Ov,q - \/30i,,q )i^sg{p )syCzCx - CyS^Sx) 

+Oxy{p )Ou,q {}Sg{p )SzCxCy - CzSxSy)) 

+ C2,i^{{Oyz{p)Oyz,q + Ozx{p)Ozx,q + Oxy{p)Oxy,q){cxCyCz + isg{p)sxSySz)) 
+ C2,b-^{pyz{p ) {Ozx,q i}sg{p )szCxCy - CzSxSy) + Oxy,q iisg{p )syCzCx - CyS^Sx)) 

+(cyclic permutation of x, y and z)^ . 

We can derive the quadrupolar interaction between the nuclear quadrupole of ^^O and the 
triple q AFQ of Np 5f electrons as 

= C2,2{-isg{p )) ( - Ouip ) + V30,{p ))Oyze'Q 



(2) 



^^qq(P , C2,2) 



+ Ou{p )OxyG 



iQ ,p 



Ou{p)-V30^{p))Ozxe 



(3) 



where the quantity C2,2 is proportional to C2,2 and the order parameter of AFQ Ops) or 
square of the order parameter of AFO (T^)^, i.e. 6*2,2 C2,20r5 C2^2{T^)^. This interaction 
originates from the coupling caused by the T^^ symmetry, such as s^CxCy term in cq.(2). At 
the p I = i(3a,3a,3a) site, the factor —isg{p i)e'Q •^P 1 takes 1 for all v. Therefore it gives 
the following expression 



Hqaip 1,02,2) 



-G 



2,2 



Ou{pl){-^Oyz+Oxy 



':Ozx) 



+ ^0,{p,){Oyz 



Oz 



(4) 



If Oyz = Ozx = Oxy = holds, this interaction is canceled out on the p i site. 

On the other hand, at the site p 2 = ^(So, 3a, a), the factor take —1,-1 and 1, respectively, 
for the terms characterized by the vectors, Q y and Q ^. Then we have 



^qq(P2,C2,2) 



G 



2,2 



G 



Oyz - O: 



1 



xy 



:0z 



:Ov{P2){Oyz - Ozx) 



= C2,20„(p2)(-20r5 



(5) 



Therefore the quadrupolar field remains on this O atom. A similar calculation can be carried 
out for other sites and the result is shown in Table III. (The sign of —\sg{p )eQ is given 
in the Table II.) We note 0„ of the ^''O nucleus (/ = 5/2) is proportional to 3/f — /(/ + 1), 
and -\{0u + \/30„) oc 3/| - /(/ + 1), and -i(0„ - \/30^,) oc 3/^ - /(/ + 1). The triple q 
AFQ order causes 4 different quadrupolar fields on the eight O sites in the fee cube: a pair 
of cubic sites (with zero quadrupolar field) and 3 pairs of uniaxial-symmetry sites ( each of 
which has principal axis along one of the x,y and z axes). The ^^O nucleus on the non-cubic 
site will show the the quadrupolar splitting of the spectrum. The appearance of the cubic and 
non-cubic site in the triple q AFO state can be understood in an intuitive way as explained 
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Table I. Product of octupolar moment arranged as the quadrupolar moments. 

5 — ^zx — J-z J-x ^xy — J-x J-y 

-p /n 1 f_T'f^Tl^ rpl3rpl3 I r,rpl3rpfi\ ^ 1 / rpP rp0 rpl3rpf3\ 

^ 3 '^u — V -Lx -i-x -i-y -'-y ~r ■^-L z z ) — ■'■x -'-y -'-y ) 



in Figs. 1 and 2. The colored cubes are the cubic sites and their neighboring three sites have 
principal axis along the x, y and z axes. 

The quadrupolar field of the type Oyz, Ozx and Oxy on O atoms is not induced in the triple 
q ordering. 




Fig. 1. Schematics for the triple q ordering of AFO. The moments on Np ions are denoted 
by arrows. The O sites are located at the centers of colored and colorless cubes. The unit cell of 
the ordered phase is composed of four Np and eight O sites. See also Fig. 2 for details. 



3. Magnetic Field Effect in Np02 

The applied magnetic field induces AFM moments with the same wave vector. They can be 
derived by using the Table III of ref. 11: 

Jx{R) = '-^^^{OxyHye^Q^R +OzxH.^Qy^) 

= ''-^^^Or.^Hy.^Q^R +Hz^QyR), (6) 

Jy{R) = '-^^Or^Hx^Q^R +Hz^Q^R), (7) 

Jz{R) = '-^^Or^iHx^Q.^ +Hy^Q^R), (8) 
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Fig. 2. Environments around O sites located at the center of cubes. There are eight nonequivalent O 
sites, (a) Two of them are subject to cubic fields. Four arrows point the center of the blue cube, 
whereas they go out from the center of the red cube. Both keep four trigonal axes, so that a cubic 
symmetry is retained there, (b) Two O sites in the colorless cubes are neighboring with blue and 
red colored cubes along z axis. Thus, the trigonal symmetry is broken there, and they are subject 
to tetragonal fields. There are other four O sites neighboring with colored cubes along x and y 
axes, respectively. These can be represented by the corresponding rotations of the figure (b). 



Table II. The sign on eight O sites and the hyperfine field on O sites. The column 

of Q J, shows —isg{p )e^^ . The last column shows i?hf , eq.(9) except for the factor Ci^2- 



site(p ) 




Qy 


Qz 


total 


§(3,3,3) 


1 


1 


1 


21 H 


§(3,3,1) 


-1 


-1 


1 


—2IzHz 


§(1,1,3) 


-1 


-1 


1 


—2IzHz 


§(1,1,1) 


1 


1 


1 


21 H 


§(3,1,3) 


-1 


1 


-1 


2IyHy 


§(3,1,1) 


1 


-1 


-1 


—2IxHx 


§(1,3,3) 


1 


-1 


-1 


2IxHx 


1(1,3,1) 


-1 


1 


-1 


— 2IyHy 



where h{T4^,T^) is a proportionality constant. Substituting this expression into eq. (1), we 
obtain the hyperfine field as 



^^hf(p,ci,2) = Ci,2is5(p) 

+IyHy (e'Q + e^Q ) + IzHz (e'Q + e'^ yP 



(9) 
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Table III. Quadrupolar field on eight O sites. Ou and Oy arc the quadrupolar operators of O nucleus. 
The last column shows the result for Oy^ = O^x = O^y = Ops 



site(p ) Ou Oy total 



2(3,3,3) 


In 


- O -x y 


+ 


In 

2OZX 


2 


'^yz 


1 /n ^ 

+ Ozx) 


U 


§(3,3,1) 




Oxy 




2 ^zx 


2 


(Oyz 


-Ozx) 


-20„Or5 


§(1,1,3) 


2^y^ 


Oxy 




kOzx 


2 


(Oyz 


-Ozx) 


-20„Or6 


§(1,1,1) 


2^yz 


Oxy 


+ 


2 ^zx 


V3( 
2 


-Oyz 


+ Ozx) 





§(3,1,3) 


2^yz 


~l~ Oxy 


+ 


2 "-^2:0; 


2 


(Oyz 


+ Ozx) 


{Ou + \/30,)Or5 


§(3,1,1) 


2^yz 


+ Oxy 




hO^x 


V3( 
2 >^ 


-Oyz 


Ozx) 


(0„ - ^/SO.Ors 


§(1,3,3) 


2^yz 


~l~ Oxy 




hOzx 


2 V 


~Oyz 


Ozx:) 


(0„ - v^O„)Or5 


§(1,3,1) 




~^ Oxy 


+ 


2'-^2a; 


2 V 


-Oyz 


Ozx) 


(0„ + V30^)Or5 



where the quantity Ci^2 is proportional to ci^2 and the AFQ order parameter, i.e. Ci^2 ^ 
ci,20rsb(T4^, Fs). This type of interaction appears through the coupling term with the constant 
ci^2 in eq. (1). If one considers only the term of the type I ■ J, which has the coupling constant 
ci^i in eq. (1), the interaction expressed by eq. (9) will not appear. In addition, this interaction 
originates from the coupling caused by the T^ symmetry, such as SzCxCy term in eq.(l). 

At the Pi = ^(3a, 3a, 3a) site, the factor —\sg{pi)e^^'^Pi takes 1 for all v. Therefore 
we have Hii{{p 1,01^2) = — 2Ci^2-^ • H. On the other hand, at the site P2 = ^(3a, 3a, a), the 
factor takes —1,-1 and 1, respectively, for the terms having Q y and Q ^. Then we have 
-f^hf (p 2, '^1,2) = 2Ci^2-^z-ffz- Clearly the hyperfine field on the pi site is different from that 
on the p 2 site. We find in this way that the eight O sites split into four inequivalent sites 
also in hyperfine interaction, as summarized in Table II. When the magnetic field is given by 
{Hx, Hy, Hz), four types of the hyperfine field are induced on the O atoms: 2Ci^2{Hx, Hy, Hz), 
2Ci,2{0,0,-Hz), 2Ci,2{0,-Hy,0), and 2Ci,2(-^x, 0, 0). 

For polycrystallinc samples wc have to take an average of the splitting over the direction. 
Then, wc find that the first group docs not show any broadening of the resonance, because 
the hyperfine field is always induced in the direction of the applied field. On the other hand, 
the hyperfine field usually deviates from the direction of the field in the other three groups. 
This will cause the broadening of the resonance spectrum with the line shape of the powder 
pattern of the axial symmetry crystal. 

Let us remember that the magnetic field usually induces the AFO moments of the types, 
Tl^ ( whose wave number vector is not equal to that of the original one) and T^yz, as shown in 
Table III of ref. 11. The hyperfine field on the O atoms is induced also through these moments 
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as can be seen in eq. (1). They give the hyperfine field of the type: 
^hf (P , ci,3) = -Ci,3is5(p ) [hH, (e'Q - e'Q yP ) 

+IyHy (e'Q ^P - e'Q ^P ) + hH, (e'^ yP -e'^^P)], (10) 

Hm{p,c^a) = Ci,4isg{p)[e'Q ^P la^H^ + e'Q yP lyHy + e'Q ^P I,H,y (11) 

where Ci,3 oc ci^aOps and Ci^ oc ci^Or^- The former gives the hyperfine field — 2Ci,3(0, 0, 0), 
—2Ci^3{Hx,Hy,0), —2Ci^3{0,Hy,Hz), and —2Ci^s{Hx,0, Hz), whereas the latter gives 
Ci,4{Hx,Hy,Hz), CiA{—Hx,—Hy,Hz), Cia{Hx, —Hy, —Hz), and CiA{—Hx,Hy,—Hz). 

To apply this theory to Np02, we have to sum up these three contributions. Even after the 
summation we still have a sharp line and a broad line; the shift of the sharp line and the field 
induced part of the width of the broad line have comparable magnitude; they are proportional 
to the strength of the magnetic field in low field region. 

Our theory has already predicted that the quadrupolar splitting with uniaxial symme- 
try appears on the non-cubic O site in the ordered phase of Np02. Even at zero field, 
this quadrupolar splitting exits. By the application of the magnetic field in the AFO state, 
the AFQ is also induced as given in the table IV. We can see that the induced quadrupo- 
lar moments of the T5 type: Oyz{R ) = ^^^{HzT^e'^ - HyT^&Q ), Ozx{R) = 
^-^{H^Tz^^Q zR -HzTU^^R) ^nd Oxy{R) = ^-^{HyT^e^Q -H^Ti^Qy^) 
cause the quadrupolar splitting. Here c(r5,r5) is the proportionality coefficient in Table IV. 
Substituting these results into eq. (2), we obtain, 

H^^{p , C2,5) = C2,5is5(p ) [Oyzip )H, [e^^ ^P - e'^ yP ) 

+OzAp )Hy (e'Q ^P - e'Q ^P ) + 0,,(p )Hz (e'^ yP - e'^ ^P )] (12) 

where €2,5 oc C2,5c(T5, r5)T''. The phase factor of this expression has the same form as that of 
eq. (10); thus it will be canceled out on the cubic sites. On the non-cubic sites, the quadrupolar 
interaction terms of the type C2,5[OyzHx + OzxHy], C2,b[OzxHy + OxyHz] and C2,5[OxyHz + 
OyzHx] are expected. 

The broadening due to the quadrupolar splitting is caused at non-cubic sites when the 
average is taken over the direction for polycrystalline samples. Even at zero field, the broad- 
ening due to the quadrupolar splitting will exist. Thus, in the low magnetic field region, the 
broadening of the spectrum originating from the non-cubic site will be the sum of a constant 
term and the term proportional to a strength of the magnetic field. 

The field-induced AFO/AFM moments are proportional to both the magnetic field and the 
AFQ moment. Therefore they are proportional to the square of the primary AFO moment. 
Another type of broadening due to the field induced hyperfine splitting is also expected. The 
AFM and the AFO moments, which are directly proportional to the primary AFO moment 
and the square of the magnetic field can be also induced. For example, from the T5 — 2 column 
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Table IV. Quadrupolar moment as a product of magnetic field and octupolar moment. 



Quadrupole moment Fs r2 



Oyz 




j={-HyTz^ + HzTi) 


HxTxyz 


Ozx 




-^{-HzTi + HxT,^) 


HyTxyz 


Oxy 




j^{-HxTi + HyT^) 


HzTxyz 


Ou 












-HxTi - HyTy^ + 2HzTz^) 





Table V. r4 type moment as a product of magnetic field and octupolar moment. 



r4 moment 




- 1 




- 2 




rpa 

y 

rpa 

z 




- H^z)Ti 
-Hl)T^ 


Hx{HyTy 
Hy{HzTz 

Hz{HxTx 


- HzTz^) 

— HxT'x) 

-HyTi) 


'\/2,Hy HzTxyz 
\/^Hz Hx Txyz 
y/^Hx Hy Txyz 



of Table V and eq. (1), we have the interaction 

Hm{p ) = Ci,2isgip ) [ixHyHz (e'Q yP + e'<3 .P ) 

+IyHzHx (e'Q + e'Q ^P ) + hHxHy (e'Q ^P + e'^ )] , (13) 

where Ci^2 ci^2T^. This term has the largest contribution to the cubic site, and the splitting 
due to this term depends on the field direction as H'^hxhyhz- After averaging over the direction, 
this effect causes a broadening of spectrum of 0{H^). However, it still gives two contributions 
with the intensity ratio 1: 3. The equations (9) ~ (11) must contain the 0{H^) terms coming 
from the change of the order parameters. 

Finally we conclude that the shift of the sharp line, Ug depend on the field strength as 
LOg = ai^^i^ + 0{H^). The width of the sharp line, A^. will be given as a higher order function 
of i^, i.e. As = 0{H^). The width of the broad line, A5, is given as A;, = c^"^ + 0^"^ H +0{H'^), 
while the shift of the broad line, cjf, is loi, = a\^^ H + 0{H'^). The constant c^"''' is caused by the 
quadrupolar splitting. If the hypcrfinc interaction is the main mechanism, the constants 0,5 , 
a^^\ will have same order of magnitude, a'^^ and a[^'' are related to the angle average of 
the axial symmetric sites. These four constants arc proportional to the AFQ moment; thus 
they must have almost the same temperature dependence. If the magnetic field induced AFQ 
in the AFO state, eq (12), has appreciable contribution, its effect will appear mainly in the 
Ol'"^ term. This contribution is directly proportional to the AFO order parameter. 
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T5 type moment as a product of magnetic field and octupolar moment. 



Fs moment 




- 1 




r5-2 


r5-3 




-U^^i - 


Hi 




-H^{HyT^ + H,tS) 






-U^^l - 


Hi 




-Hy{H,T^ + H^Ti) 






-u^^i - 


Hi 


- Hl)Ti 


—Hz{HxTi + HyT^) 





4. Hyperfine Interaction with Multipolar Moments in CeBg and Cea.Lai_a.B6 

CeBg has the crystal structure of CaBg type. The Ce ions form the s.c. lattice. Let us denote 
the B ion site (|,|,ib'ua) as p z,± = (OjOj^^)) where n is a parameter for the position of 
B atom and a is the lattice constant of the s.c. lattice. If we follow the method described in 
Appendix C, the interaction of the nuclear spin J of B atom with the multipolar moments of 
Ce ions is given by 

i?hf(p.,±) = ^e'^^V9^^.±[4(p,,±){ci,iJ,(g)c,c, 

+i«5(P z,±y^ {Mq )sxCy + Jy{q )c^Sy^ + isg{p (q )s^Cy - T^{q )c^Sy^ 

~Cl^4Txyz{Q )SxSy"^ 

+ Ix{p ^^±)|c2,l Jx(q )CxCy + isg{p ^^±)c2,2Jz{Q )SxCy - C2,2,Jy{q )SxSy 

+C2AT^{q )cxCy + C2,5T^{q )sxSy + [sg{p ^^±)c2,6Tf (q )sxCy + [sg{p ^^^)c2,7Txyz{q )cxSy^ 

+ lyiP z,±){c2,lJy{q )CxCy + isg{p ^^±)C2,2 )SyCx - C2,2,Jx{q )SxSy 

+C2,4Ty{q )cxCy - C2,5T^{q )sxSy + isg{p ^^^)c2,%T^{q )syCx 

+^sg{p ^^^)c2,7Txyz{q )c^s^}] ■ (14) 

Here c^j's are coupling constants, and sg{p ^ j.) = itL The quantity Ci,{si,) represents 
cos{qi,/2) (sin(gi,/2)). The result (14) agrees with the result given in ref. 9. 

The hyperfine interaction of the B nuclei at the (ibua, |, |) site, which is denoted hence- 
forth as p x^± = (±^! 0, 0), is obtained by 90° degree rotation around the y-axis: {Ix, ly, Iz) — *• 
[—Iz,Iy,Ix), {Tx-i~Ty,Tz) > (Tz , —Ty , —Tx) , Txyz ^ ~Txyz- At the same time the 
wave vector is changed as {qx,Qy,Qz) {—Qz^Qy^Qx) except for the factor in the func- 
tion e 2 Q'iHx,±_ The dipole operator [Jx,Jy,Jz) follows the same transformation 
to that of the operator {Ix,Iy,Iz)- Similarly, the hyperfine interaction for the B on the 
(|,ibua, |) site, which is denoted as p y^± = (0, ±^,0), is obtained by a rotation around 

the X-axis: {Ix, ly, Iz) ^ {,Ixi~Izi Iy)i {'^x :~Ty ,Tz) > {—Tx ,Tz i~'^y)i '^xyz ^ ~'^xyzi 

{(ix,qy,qz) {qx,-qz,qx)- 

If we assume the Q = 7r(l, 1, 1) order as done in refs. (7) and (28), the quantities Si, and Ci, 
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as follows: 



1, and Cx = Cy 
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0. Then the hyper fine field is expressed 



Hm{p,,±) = ±[/.(p,,±)(-Ci,4r,,,(Q)) 

+Ix{p ,,±) ( - C2,3Jy{Q ) + C2,^Ty^{Q 

+Iy{p ( - C2MQ ) - C2,,Ti{Q )) 
+7,(p,,±)(-Ci,iJ,(0)) 

+Ix{p (C2,l J.(0) + C2,4Ti^(0)) + Iy{p (C2,l J,(0) + C2,4T^^(0))(]l5) 

Here T^{Q ), J^iQ ) and Txyz{Q ) arc the thermal average of the ordering parameters with 
wave vector Q . The quantities for Q =0 are the uniform components. The factor Cj^- is a 
constant proportional to Cij. Even in the uniform term wc have the hyperfine interaction due 
to the octupolar moment T^ifS). Such term will be induced by the application of the magnetic 
field under the existence of the uniform component of quadrupolar moment. It is given by 
a different combination of the quadrupolar moment and the magnetic field from that of J^^q 
term.^^ Hereafter we neglect the Q = terms. The hyperfine interaction for p ^ _|_ and p y^^ 
site is given , 

i7hf(p.,±) = ±[lxiPx,±)i-CiATxy,{Q)) 

+h{Px,±){ - C2MQ ) - C2,5Ty^iQ )} 

+Iy{P x,±){ - C2,3MQ ) + C2,,T^^iQ )}] , (16) 



Hhf{p y, 



±) 



-- ±[lyiPy,±)i-CiAT.yziQ)) 

+ix{p - C2,3MQ ) - C2,5rf (Q )} 

+4(P y,±){ - C2,3Jx{Q ) + C2,,T^{Q )}] . (17) 

If the ordering wave vector q satisfies cos{qy/2) 7^ cos(g^;/2), a mixing of Jx{q ) with T^^q ) 
through the nearest site interaction on the s.c. lattice^^'^^ will be caused. In addition, the 
ordering with sin(giy) 7^ causes mixing of dipole and type moment through the next 
nearest site interaction.^-*^ The ordering of the wave vector with Q = 7r(l, 1, 1) seems to be 
a plausible candidate when one considers the pure spontaneous AFO. The n.n. interaction 
of Fs moment will have comparable magnitude to that of the quadrupolar and magnetic 
interaction.^^ 
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5. Dependence of Hyperfine Splitting on the Field Direction in Cei—jcLajcBe 

Let us assume that the ordering with Q = 7r(l, 1, 1) occurs. Then, the hyperfine field 
splitting A(2;) of B-ion pair at (|, |,±?/a) is given as 

A(^) = -2Ci,4h,T^y,{Q ) + 2C2,5{h^Ty^{Q ) - hyT^{Q )) 

-2C2,3{hxJy{Q ) + hyUQ )), (18) 

where {hx,hy,hz) denotes the direction of the magnetic field. Henceforth we consider the 
ordering of three types: 

(i) ri^(Q ) = Ty^iQ ) = (Q ) = ^rfi,, 

(ii) TiiQ ) = -TiiQ ) ^ ^rfj„,rf (Q ) = 0, 

(iii) r.^(Q ) = Ty^iQ ) = -^Tfio, and Tf (Q ) = ^rf^^- 

In the case of TiiQ ) = T^yiQ ) ((i) and (iii)), the expression (18) depends on the direction 
of the magnetic field as sin^sin((/) — j). Here the direction of the magnetic field is expressed 
as (sin^cos^, sin^sin(^, COS0). A{z) vanishes at tp = 7r/4, and has linear dependence, as 
A(z) (X cj) — J, in the vicinity of (/> = j. In the case of Tx{Q ) = —Ty{Q ) (case (ii)), A{z) 
depends on the direction as sm9s'm{(p + j), showing the maximum at (p = 7r/4. 

When the magnetic field is applied, the order parameters are modified. The lowest-order 
effect is given by the product of the square of the magnetic field and the linear term of the 
order parameters, because they are the lowest order expression with non- uniform and time 
reversal odd in the case of AFO ordering. In Table V, we show the induced r4 AFO moment 
as a product of the Fs AFO moment and the magnetic field. This will also be proportional to 
the dipolc moment. Substituting this into the third term of eq. (18), we have 



hxJyiQ ) +hyJx{Q ) 



ai(F4,r5) 2 , , , , , 2 



= {hx - hy) [ ^''^^"" {hi + hxhy +h'y- hi) - a2{VA,Vr,)hxhy\ H^T^ , (19) 

where ai(F4,F5) and a2(r4,r5) are proportionality constants. Here we introduced Ti{Q ) = 
Ty{Q ) = . The terms proportional to the order parameter Tz{Q ) do not appear even 
when it exits. Equation (19) is also proportional to hx — hy, thus the linear dependence on 
(p — tt/A is not changed. The induced moments of Fs type can be derived also by using the 
Table VI, but the conclusion is unchanged. 

For the order of the type Jx{Q ) = —Jy{Q), the same dependence as that of Tx{Q ) = 
Ty{Q ) is expected as seen from eq. (18). The field dependence is also similar as one can check 
easily by using the results of Tables VII and VIII. However, complementary measurements of 
the splitting for pairs on the other axis, and/or detailed calculations of the direction depen- 
dence based on the microscopic model will distinguish them. If the Txyz{Q ) order is realized, 
the dependence of A(z) on the field direction is unique. This has been already noted in ref. 9. 
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Table VII. r4 type moment as a product of magnetic field and dipole moment 



r4 moment 


r4 


- 1 






r4-3 


-'■X 

rpa 

y 

rpa 

z 




H.I 
Hi 
Hi 


- Hl)Jx 

- Hl)Jy 

- Hi) J, 


— Hx{HyJy + HzJz) 

—Hy{HzJz + HxJx) 

— Hz{HxJx + Hyjy) 


T^H^-^^ 



Table VIII. r5 type moment as a product of magnetic field and dipole moment. 



moment 




1 




r4- 


- 2 


rpP 

^y 

rp(3 
J-Z 


T^^Hl - 

^^Hl- 


Hi] 
Hi. 

Hy, 


)Jx 

)Jy 

)Jz 


Hxi^HyJy 
Hy{HzJz 

Hz{HxJx 


— HzJz) 
HxJx) 

— HyJy) 



Table IX. r2(ri) type moment as a product of magnetic field and the r4(r5) moment 



r4 




Txyz — ~ sj^iHyH zJx + HzHxJy + HxHyJz) Ai = ^ 


J '^{HyHzTi + HzHxTy + HyHzTz) 



The splitting of the NMR of the B ion pair at (±ua, §, §), A(a;), is given as 
A(x) = -2Ci,4hxT^iQ ) + 2C2,5 (hzT^iQ ) - hyT.^iQ )) 

-2C2,3 (hzJyiQ ) + hyUQ )) . (20) 

The splitting of the pair at (|,ibna, |), A(y), is similarly given by using eq. (17). When the 
octupolar ordering follows Tz{Q) = Ty{Q ) (type (i)), this expression has a node when hz = 
hy and linear dependence around the direction. The ordering of the type Jz{Q ) = —Jy{Q ) will 
also give the same direction dependence. However, they can be distinguished when referring 
to the experimental result of A{z). The ordering of the type Jx{Q ) = Jz{Q ) = —Jy{Q ) does 
not have symmetry around the [111] direction. Such ordering usually does not appear. But 
the ordering of the type (i) can appear. Moreover, the splitting A{y) has different direction 
dependence between them. 

In the ordering of type (ii), the splitting A(x) has the largest value for the field in the z 
direction. In the ordering of type (iii), the splitting A(x) becomes zero for hz/hy = —2, and 
takes a maximum at hy/hz = 2. 
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6. Summary 

Wc have presented a new approach to derive the hyperfine interaction between the nuclear 
spin of Ugand atom and multipolar moments of magnetic ions. 

The most important part of this paper is the analysis on -"^^O NMR in Np02, which has been 
reported recently by Tokunaga et al. et. al}^ We have studied the ^'^O NMR spectrum for the 
triple q order of primary AFO and secondary AFQ, which was proposed by Paixao et. al.^ At 
zero field the triple q AFO does not produce any hyperfine field on the O sites. However, the 
secondary triple q AFQ order causes 4 different quadrupolar fields on the eight O sites in the 
fee cube: a pair of cubic sites (with zero quadrupolar field) and 3 pairs of uniaxial-symmetry 
sites (the principal axis: x,y or z axis). When the magnetic field is applied, the AFM/AFO 
is induced in cooperation with the preexisting triple q AFQ order. As a result, the NMR 
spectrum from the cubic site gives a sharp line with a shift proportional to the strength of 
the magnetic field because the hyperfine field is isotropic. On the other hand, the spectrum 
from three kinds of non-cubic sites gives the same shape for polycrystalline samples. Namely, 
it has a broadening due to the quadrupolar splitting and the magnetic field induced part 
caused by the anisotropic hyperfine field. The experimental observation by Tokunaga et al. is 
summarized as follows: 

(I) splitting to sharp and broad lines with the ratio 1:3 in the ordered phase, 

(II) the magnetic field dependence of the shift of the sharp line, 

(III) the magnetic field dependence of the width of the broad line, 

(IV) relation between the shift and the magnetic field induced part of the width. 

Those features are consistently explained by our scenario. The present calculation and the 
experimental result in 19 strongly support to the triple q ordering model of Np02. 

We have found theoretically that in the AFO state there are unique coupling terms such 
as the one causing the quadrupolar splitting proportional to the magnetic field and/or the 
hyperfine field splitting proportional to the square of the magnetic field. The NMR experi- 
ment for single crystals is highly desirable. On the theoretical side a microscopic calculation 
including the effect of finite magnetic field is desirable. In addition, although our phenomeno- 
logical approach is certainly useful, a microscopic calculation of the transferred hyperfine 
interaction based on a microscopic model^^ is also desired. Such microscopic calculations may 
modify quantitatively the results of the present paper, but we believe they will not change 
our conclusion drastically. 

In the second part we have derived the hyperfine interaction between the B nuclear spin and 
the multipolar moments of Ce ion in CeBg and Cei-x-Lax-Bp. It has been applied to discuss 
a possibility to identify the octupolar order parameter, which was proposed for phase IV 
of Cei_a;La^.Bg. We assumed the pure octupolar order of the r5(r'^) type with wave vector 
Q = 7r(l, 1, 1) under uniaxial stress of the [1, 1, 1] direction. If it is realized in Cei-^jLaajBe, 
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the hyperfine field splitting of B nuclear spin should show characteristic features dependent on 
the type of order. For example, if the order of type Tx = Ty = Tz occurs, the splitting l^{z) 
of NMR vanishes when the direction of the magnetic field crosses the (110) plane. Therefore, 
in principle, the octupolar ordering can be identified from the direction dependence of the 
hyperfine splitting. Magishi et. al. have made an NMR experiment on the phase IV. ^'^ Though 
the observed lines are sharp enough in phase I (normal phase) , they overlap and become 
broad in phase IV. They speculated from the experiment that an AF magnetic ordering of 
incommensurate wave vector is realized. The NMR experiment under a strong uniaxial stress 
along (1,1,1) is desirable. Although the pure AFO ordering has been assumed for phase 
IV, Cei_a;Laa;Bg is actually a disordered system so that the AFO is accompanied inevitably 
by AFM clusters because of low local symmetry. This effect should show up in the neutron 
scattering experiment. In this connection we note that strange insensibility of Ce ions to the 
magnetic ordering in Ndi_a;Cea;B6 may be related to this problem. ^^'^^ 
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Appendix A: Derivation of Invariant Coupling Form of Hyperfine interaction in 
Np02 

Let us denote the O site ^) as p ^. Wc consider the invariant couphng form between 

the nuclear spin / of ^^O on this site and multipolar moments of f-electrons of Np ions at four 
nearest neighbor sites, (2a, 2a, 2a), (2a, a, a), (a, 2a, a) and (a, a, 2a) of the cube of the f.c.c. 
lattice. For simplicity we denote these Np ions as (111), (111), (111)) and (111), respectively. 
These ions have the symmetry around the O atom. In a previous paper, the invariant 
form is derived by considering the pair interaction and rotating the pair by the symmetry 
operation.^ This method certainly gives the correct interaction form, but as noted in ref. 11, 
the calculation becomes relatively easier if one uses the symmetrized combination of multipolar 
operators. The nuclear dipole operators Ix, ly and form the 3-dimensional Ti representation 
of Tj.^^ Following the method to make the symmetrized molecular orbital, we can construct 
combinations of and Fs and F2 type odd power operators of the Np ions, which have Ti 
representation. 

Then, the invariant form of the hyperfine interaction of bilinear type is given as follows: 
Hhi{p 1) = IxiP 1) [ci,l;^|(Ja;)(lll) + ('4)(lll) + (■4)(IlI) + ('4)(111)} 

+Cl,2^|(-/y + Jz)){Ul) + {-Jy - Jz)(in) + {-Jy + <^z)(ill) + (<4 - <^z)(lll)} 

+ci,3^{{TyP - ))(iii) + {-Ty^ + )(iii) + (-7? - )(ni) + (T^ + T^,, 

+Ci^4:-^^{Txyz){Ul) + {Txyz)(lll) + {—Txyz)(lll) + {—Txyz) (111)^ 

+(cyclic permutation of x,y and z) . (A 

Here Cij are the coupling constants. Next we make the Fourier transformation for multipolar 
operator of Np ions, 0-y{Ri) = SqO-y^g e'^ where R ^ denotes the site of Np ion, 
and 7 denotes the symmetry of the multipolar operator. 

The interaction for nuclear spins at (|,|,^), (^,|,|) and (f,^,§) can be obtained 
by translating the site by (— a, — a, 0), (0, — a, — a) and (— a,0, ~a), respectively. While the 
interaction of nuclear spins at (^, ^, f ), (f , |, f ), (^, f , ^) and (|, ^, ^) can be obtained 
by 90 degree rotation around z axis of the T^ and subsequent translation of the site by 
(0,0, —a), (—a, —a, —a), (0, — a, 0) and (—a, 0,0), respectively. By the rotation, the operators 
are changed as follows: (4, /y, h) {ly, -Ix, h), {Jx, Jy, Jz) ^ {Jy, -Jx, Jz), {t!^ ,T'^ , Tf ) 
{-T^,T^, -Tf ) and T^y^ ^ -T^y^. The sites are changed as (111) ^ (111), (111) ^ (111), 
(111) ^ (111) and (111) ^ (111). 

The Fourier transformation is carried out in the same way as was done for eq. (A-1) of the 
p I site. Note that the same expression can be obtained by transforming the wave vector in 
the Fourier transformed equation for p i. e. {qx, Qy, Qz) {Qy, —Qx, Qz) except for the factor, 
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Q^Q P . Of course, the operators are changed in the way noted above. 

Appendix B: Quadrupolar coupling in Np02 

Let us consider the quadrupolar interaction between the r3-type nuclear quadrupolar mo- 
ment of O nucleus and the Fs-type quadrupolar moment of Np ions, 

Hc^dip 1,C2,2) 

= C2,2 Ou{p i)--^^-^^{Oyz + Ozx - 20a;y)(iii) + (Oyz - Ozx + '^Oxy){iii) 
+{—Oyz + Ozx + 20xy)(iii) + {—Oyz — Ozx — 20a;y)(ni)| 

+ Ou{P l) Ozx)(\l\) + {-Oyz - Ozx){lll) 

+ {Oyz+Ozx){lll) + {Oyz-Ozx){lll)^ ■ (B-l) 

Carrying out the Fourier transformation, we get 

-f^qq(Pl,C2,2) 

8 r r 1 ~i 

= e'^ P 1C2,2 ^^ |-(-0„(p l) + V2,Oy{p l))Oyz,q {-iSxCyCz + CxSySz)j 

+{Ou(p )Oxy,q {-iszCxCy + CzSxSy)^ 

+ |^(-0„(p i) - VS0^{p i))Ozx,q {-^SyCzCx + CySzSx)^ . (B-2) 
This corresponds to the interaction term with the coupling term C2,2 in eq. (2). 

Appendix C: Derivation of Invariant Coupling Form in CeBe 

Let us consider the invariant coupling form between the nuclear spin I of B on (|, ^,ua) 
and multipolar moments of f-electrons of Ce ions at (a, a, a), (0, a, a), (0,0, a) and (a, 0, a) of 
the cube of the s.c. lattice. For simplicity we denote these Ce ions as (111), (111), (111), and 
(111), respectively. 

These ions have the C4v symmetry around the z axis. The nuclear spin operator Iz of B 
belongs to the A2 representation, and {Ix,Iy) form a two-dimensional E representation. 

We can construct combinations of multipolar operators of the Ce ions, which belong to A2 
and E representations. Then, the invariant form of the hyperfine interaction of bilinear type 
is obtained as follows: 

Hhf{^, ^,ua) = Iz ci,i-^^{Jz) (111) + {Jz) (111) + {Jz) (111) + {Jz) (111)} 

+Cl,2-^\^{Jx + Jy)(lll) + {-Jx + Jy)(lll) + {-Jx " Jy)(lll) + {Jx " •A/)(iii)| 
+ci,3^{(rf - Ty%n) + {-Ti - Ty%,,) + (-Tf + Ty%-,,) + {T^ + Ty%-,,)} 

+Cl,A-^^{Txyz)(lll) + {—Txyz)(lll) + {Txyz)(ni) + (" ^a:j/z)(lll) | j 
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+ I 



C2,l^|(Ja:)(lll) + (<4)(111) + (-4)(111) + (<4) 
+C2,2;^|(^z)(lll) + {-Jz){lll) + {-Jz){lll) + (•^z)(lll)} 
+C2,3-^|(-'2/)(lll) + (--'y)(lll) + {Jy){lll) + i-Jy){lll)\ 

+C2,4^{(T;^)(iii) + (7f )(iii) + (Tf )(iii) + (Tf )(ni)} 
+C2,5^{(-r,^)(iii) + (r,^)(in) + + {T^\ili)] 

+ + + (rf )(iii)} 

+C2,7-^^{Txyz){lll) + {Txyz)(ll\) + (-ra;y2)(iii) + (-ra;j,2)(iii)|j 
+ ly C2,l-^^{Jy){lll) + (Jj/)(111) + (-'2/)(lIl) + ('^2/)(lll)} 

+C2,2;^|(^2)(111) + (-'z)(lll) + {-Jz)(lll) + (--'2)(111)} 
+C2,3;^|(-^3;)(111) + + + (-•4)(111)} 

+C2,4^{(r,^)(in) + (r^)(in) + (T?)(in) + 
+C2,5^{(7;^)(iii) + )(iii) + + )(in)} 

+C2,6^{(-ri^)(iii) + (-Tf )(in) + (Tf )(IIi) + (7f )(iii)} 

+C2,7;^|(ra,y2)(iii) + + {-Txyz)(lll) + {Txyz)(lll)^ ■ 

Here Cj y are coupling constants. By carrying out the Fourier transformation, we get the 
interaction form given in eq.(14). 

The interaction for the nuclear spin of B on (|,|,— «a) can be obtained by a mirror 
operation which moves B on (|,|,Ma) to By this transformation the mul- 

tipolar operators are changed as follows: {Ix,Iy,Iz) i—Ix,—Iy,Iz), iJx,Jy,Jz) 
{-Jx,-Jy,Jz), {Ti,T^,T^) {-Ti,-T^,T^) and Txyz Txyz- At the same time, the 
sites are changed as (111) ^ (111), (111) ^ (III), (111)^(111), and (111) -> (111). 
The Fourier transformation is carried out in the same way as was done to derive eq. (14). 
Note here that the same expression can be obtained by transforming the wave vector in eq. 
(14) as {qx, Qy, Qz) [Qx, ly, —Qz) except for the factor, e 2 e'^ ^ . 



(C-l) 
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